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Calorimetric Determination of the Half -Life of Polonium 

D. C. Ginnings, Anne F. Ball, and D. T. Vier » 

The heats of radioactivity of four samples of polonium have been measured with a 
Bunsen ice calorimeter over a period of about seven months. With samples ranging in 
initial powers from 0.17 to 1.4 watts, the half -life values calculated from these measurements 
were found to agree within 0.1 percent, or the equivalent of 0.0003 watt, whichever was the 
larger. The results with the sample with the largest power gave a half -life value of 138.39 
days, with an uncertainty of 0.1 percent (0.14 day). This is in agreement with the value 
of 138.3 days (±0.1 percent) reported by Beamer and Easton, who used a different calorimet- 
ric method. 



1. Introduction 

Of all published determinations of the half-life of 
polonium, probably the most accurate is that hy 
Beamer and Easton [l], 2 who used a calorimetric 
method and observed a sample for as long as 97 days. 
Their method consisted in measuring the temperature 
difference between a container with the sample of 
polonium and an identical container that was empty. 
Both containers were surrounded by an isothermal 
jacket, and the relation between temperature differ- 
ence and power was determined by means of electric 
calibration experiments. The value of the half-life 
was found to be 138.3 days (±0.1%). A Bunsen 
ice calorimeter is also suitable for measuring radio- 
active power (and half-life) of radioactive materials. 
Because this method is entirely different from that 
used by Beamer and Easton, and because much larger 
samples of polonium were available, it seemed 
worth while to measure the half-life of several sam- 
ples that were suitable for measurements with the 
ice calorimeter. 

The Bunsen ice calorimeter seems ideally suited for 
several reasons to the measurement of the heat 
evolved by curie or multicurie quantities of alpha, 
weak beta, and with modification [2], other radio- 
active materials. First, the heat leak is small (per- 
haps 0.0002 w). Second, the calorimeter requires 
very little attention during the measurements, 
usually a few minutes every hour or two for ice-bath 
replenishment. Third, no electrical or temperature 
measuring instruments are required because the 
measurement of heat requires only weighing of mer- 
cury. Fourth, the calibration factor of the ice 
calorimeter is a fundamental physical constant, 
which has been determined to about 0.01 percent by 
electrical calibration experiments [3]. 

On the other hand, the ice calorimeter in its present 
application has an inherent variation, or error, which 
is comparable with its heat leak. While this does 
not limit the precision of experiments involving about 
1-w power, it would limit the precision of measure- 
ments of much smaller powers. Other ice calorim- 
eters designed specifically for small radioactive 
powers [2] have proved sensitive to as little as 
0.00003 w. An improved ice calorimeter has been 

1 Los Alamos Scientific Laboratory, University of California. 

2 Figures in brackets indicate the literature references at the end of this paper. 



in use at the National Bureau of Standards for some 
time in the measurement of heat capacities at high 
temperatures [3]. Although this calorimeter was not 
intended for measurements of radioactive power, it 
seemed suitable for measurements on certain samples 
of polonium (alpha emitter) furnished by the Los 
Alamos Scientific Laboratory, which initially de- 
veloped powers in the range of 0.17 to 1.4 w. By 
measuring the decay of radioactive power of a sample 
of polonium over time intervals comparable with its 
half-life, values were obtained of its half-life. 

2. Experimental Details 

The ice calorimeter used in these experiments is 
shown in figure 1. It is the same calorimeter pre- 
viously described [3]. An ice mantle (I) was frozen 
around the central calorimeter well (A) on a system 
of copper vanes (F) designed to increase the area of 
contact of the central well with the ice. The sample 
was suspended in the bottom of this well by a small 
wire that was pushed over against one side of the 
well by the gate G. The heat developed in the 
sample melted some of the ice, thereby decreasing 
the volume of the ice-water system and causing mer- 
cury to be forced into the calorimeter from beaker 
B. A small flow upward of dry helium in the calori- 
meter well was maintained to prevent condensation 
of water vapor from the room in the calorimeter and 
also to increase thermal contact between the sample 
and the calorimeter. 

The polonium samples were sealed in glass or 
metal containers, which were enclosed in brass outer 
containers made to fit closely to the central well of 
the ice calorimeter. Because the radiation from 
polonium is essentially all of the alpha type, the 
containers were completely effective in converting 
all the radioactive energy into heat within the con- 
tainer. A sample was lowered into the calorimeter 
(with the two platinum shields above it, as de- 
scribed elsewhere [4]), and about 10 or 15 min was 
allowed for it to come to essentially a steady tem- 
perature, the valve V being open. Then the mer- 
cury meniscus was adjusted to the upper part of 
the scale, C, by partly evacuating reservoir R. 
Valve V was immediately closed, thereby causing 
the meniscus at C to start to move downward due 
to the heat input to the calorimeter. The time wa 
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Figure 1. Schematic diagram of ice calorimeter. 

"" A, calorimeter well; B, beaker containing mercury; C, glass capillary; E, ice 
bath; F, copper vanes; G, gate; I, ice mantle; M, mercury; P, Pyrex containers; 
R, mercury reservoir; T, mercury "tempering" coil; V, needle valve; W, water; 
1, 2, 3, and 4, thermocouple junctions. 

recorded Hvhen the meniscus passed an arbitrary zero 
on this scale, the beaker of mercury (B) was re- 
placed by a weighed beaker of mercury, and valve 
V opened again. At the end of an experiment, the 
valve was closed again and the time recorded when 
the meniscus again passed its zero. The difference 
in the weights of the beaker gave the weight of mer- 
cury sucked into the calorimeter in a measured time. 
A small correction (perhaps 1 joule) was made for 
heat leak. The corrected mass of mercury was con- 
verted to energy by use of the calibration factor of 
the ice calorimeter as previously determined [3], 
270.46 abs j/g of mercury for the "ideal" calorimeter 
or 270.47 abs j/g for the actual calorimeter used. 
From the duration of the experiment (difference in 
times at which the mercury meniscus passed the 
arbitrary zero), the energy was converted into 
average power. 

The values of power calculated in the manner de- 
scribed would be correct if the sample and all other 
parts in thermal contact with the calorimeter did 
not change temperature during the experiment. 
However, as heat is put into the calorimeter, the 
thickness of the water layer between the ice and the 
central well with its vanes increases, and as there is 
flow of heat across this water layer, the temperature 
drop across this water layer must increase. As a 
result, the temperatures of the central well and 
vanes, and the sample with its container, all in- 
crease slightly during an experiment, so that part of 
the heat produced by the sample does not melt ice. 
This amount of heat depends on the power of the 
sample, the heat capacities of the various parts, and 



on the thickness of the water layer at the beginning 
of the experiment. The analysis of the correction 
for this heat is given as follows. 

Consider a system (such as the radioactive sample 
plus metal calorimeter well with vanes) separated 
from ice by a layer of water. Consider also that the 
thermal contact between sample and calorimeter 
does not change during an experiment. Let 

Q=the total amount of heat (joules) that has been 
introduced into the ice calorimeter at any 
particular time, 

P=the power (j/sec) developed by the sample, 

H=ihe heat capacity (j/deg C) of the system 
(sample plus calorimeter well with vanes), 

.4= the effective area of contact of the system 
with ice, cm 2 , 

K = the thermal conductivity of water at 0° C, 
0.0052 w cnrMeg- 1 , 

m=the average thickness of the water layer, cm, 

F=ihe heat of fusion of ice, 333 j/g. 
Then the grams of ice melted = Q/F=Q/333. 

The volume of water between the ice and the vanes 
is the same as the volume of ice melted because any 
void created by the difference in density of ice and 
water is filled almost instantaneously by water out- 
side the mantle passing through and around the 
cracks in the ice mantle. Hence, the volume of 
water layer formed = (specific volume of ice) (Q/333) 
-(1.09) (ft/333) =0.0033Q; m = 0. 0033 (QJ A); and 
the temperature difference across the water film= 
Pm/(0.0052).4 = (0.0033/0.0052) (PQ/A 2 ) = 0.63(PQ/ 
A 2 ). Thus the heat stored in the system after Q 
joules have been introduced (assuming no tempera- 
ture gradient in the metal system) = (0.63H) (PQ/A 2 ), 
and the heat stored in the water after Q joules have 
been introduced (assume linear temperature gradient 
across the film) = (4.18) (0.00330) (0.63/2) (PQ/A 2 ) = 
0.0043 (P0 2 /A 2 ), and the total heat stored=(PM 2 ) 
(0.63HCH* 0.0043Q 2 ). If we start an experiment 
with Qi joules of heat alread}^ put into the calorim- 
eter and a steady state of heat flow from the sample 
to the ice mantle, and end the experiment with Q 2 
joules of heat in the calorimeter, and the same steady 
state of heat flow, then the difference in heat stored 
= (P/A 2 )[0.63H(Q 2 -Q l )+0.0043(Qi-Ql)l 



The relative error, E, is 



E= 



~A 2 (Q 2 -Q 1 ) 



[0.63ff(Q 2 -Q 1 )+0.0043(<?!-Q?)], or 



#=— [0.6327+0.0043 (ft-f^)]. 

If E is plotted against (Q2+Q1), the slope of the 
resulting straight line is 

P 

Slope = 0.0043-^, 

and the intercept of this line with the E axis is 

HP 

Intercepts. 63 — p- 
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This analysis shows that for any given power P of 
the sample, the slope of this line is dependent only 
on the effective area of contact (A) of the metal vanes 
with the ice. This effective area depends, of course, 
on the length of the sample. The intercept of this 
line with the relative error axis depends not only 
upon the area A but also upon the heat capacity, H 
of the sample-plus- vane system. 

For any given sample with a known heat capacity 
and length, values of A and H can be calculated, 
using the dimensions and heat capacity of the cor- 
responding portion of the calorimeter well and its 
system of vanes. It is believed that a better evalu- 
ation of A and H can be obtained by electrical cali- 
bration experiments, putting in known quantities of 
electric heat into the calorimeter over a length 
equivalent to the radioactive samples to be measured. 
If the heat is distributed in the electrical experiments 
the same as in the experiments with the radioactive 
samples, the value of A obtained from the electrical 
experiments should also be the same. However, the 
value of H obtained from the electrical experiments 
differs from that for the radioactive experiments by 
the difference in heat capacities of the electric heater 
and the radioactive samples. 

The results of the electric calibration experiments 
are indicated by the circles in figure 2, which shows 
the relative error, E, plotted against (Q1 + Q2)' 
The best straight line through these circles is labeled 
"Electric Experiments" and was determined by the 
method of least squares to be £'=0.000633 + 2.992X 
10~ 8 (Q 1 + Q 2 ), where Q x and Q 2 are expressed in 
joules. The electric calibration experiments were 
made with a power of about 1.4 w to correspond 
approximately to the initial radioactive power of 
sample 4, which was the largest of the four samples 
and therefore had the largest correction for the water 
layer. The slope of the solid line (electric experi- 
ments) indicated that the effective area of contact 
(^4) between the ice and the source of heat was about 
420 cm 2 . This is roughly equivalent to a 9-cm 
length of the calorimeter well, as compared to an 
actual length of the electric heater of about 7 cm. 

While the value of A obtained from the electric 
calibration experiments agreed as well as expected 
with that estimated from the dimensions, the value 
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Figure 2. Corrections for water layer in calorimeter. 
Electrical experiments: ©, 10-11-49; ©. 10-12-49; O, 10-14-49. 



of heat capacity H (after correction for the heat 
capacity of the heater) corresponded to about a 
17-cm length of the calorimeter well, which is con- 
siderably more than expected. Some of this differ- 
ence may be due to the assumption in the theory 
that the temperature gradient in the copper vanes ls 
negligible. If the position of the line had been cal- 
culated entirely on the basis of a 9-cm length of the 
calorimeter well, it would have the same slope as the 
line given by the electric experiments but would be 
shifted to decrease E by almost 0.0003. Because it is 
believed that the electric calibration experiments give 
more reliable values of the effective A and H, the 
solid line in figure 2 was used as basis for the correc- 
tions for the experiments with the radioactive 
samples. The dashed lines shown in the figure for 
sample 4 were calculated on this basis, as well as 
corresponding lines (not shown) for the other samples. 
Values of apparent calibration factor (which is 
equivalent to 270.47 (l-\-E)) were read from the 
dashed lines for the corresponding values of (Qi + Q 2 ). 

3. Results 

The results of the measurements with the four 
samples are given in table 1. The mean time of 
experiment is given, so that all experiments in one 
group can be corrected for the radioactive decay to 
bring them all to the same time (noon, EST) on the 
reference date. The mass of mercury was corrected 
for the effect of heat leak between the calorimeter 
and its surroundings (this effect is usually only a 
few milligrams of mercury). The quantities Qi &rid 
Q 2 are as previously described. Using the quantity 
Q1 + Q2, the apparent calibration factor of the ice 
calorimeter was obtained from figure 2, from the 
curve for the appropriate sample. The product of 
this factor and the mass of mercury, divided by the 
duration of the experiment, gives the power at mean 
time of experiment. Using an approximate value of 
the half-life of polonium, this power was converted 
to power at noon (EST) on the reference date. The 
deviations from mean power are the deviations of the 
results of the individual experiments from the mean 
power of the group of experiments. The values of 
half-life given in the last column are based on the 
mean powers and reference dates, and the values 
indicated in parentheses are the authors' estimated 
uncertainties of the half-life based on both accidental 
and systematic errors. 

It will be noted that more measurements were 
made with sample 4 than with the other three sam- 
ples. This was because sample 4 had the highest 
power of the four samples, and it was believed that 
the highest relative calorimetric accuracy could be 
obtained with this sample, resulting in the best de- 
termination of the half -life of polonium. The power 
of sample 1 was so low that calorimetric errors could 
account for the low value of half -life calculated. 
The uncertainty in heat leak of the calorimeter was 
believed to be about 0.0002 w. This uncertainty, 
when compared with the power measured for sample 
1 on May 14, 1950 (0.05698 w), could explain most 
of the discrepancy with the results with the other 
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Table 1. Results of experiments 



Date of experi- 
ment 



Mean time of 

experiment 

EST 



Duration 
of experi- 
ment 



Mass of 
mercury 



0i 



02 



Apparent 

calorimeter 

factor 

270.47(1+^) 



Power at 
mean 
time of 
experi- 
ment 



Power at 
noon (EST) 



reference 
date 



Reference 
date 



Devia- 
tions 
from 
mean 

power 



Mean 
power 



Half-life 



Sample 1 



10-3-49 
10-4-49. 

5-3-50. 
5-4-50 .. 

5-14-50 



3:24 p. m. 
6:36 p. m. 
12:32 p. m. 
3:20 p. m. 

1:14 p. m. 
1:19 p. m. 
2:16 p. m. 
7:14 p. m. 



sec 





?Hg 


ffHg 


3 557.8 


2. 3048 





2 


12 367. 2 


7. 9892 


2 


10 


9 928. 


6. 3924 


10 


17 


10 126.5 


6. 5205 


17 


23 


20 354.0 


4. 5269 


24 


29 


19 440. 


4. 2913 


29 


34 


21 909. 5 


4. 6221 





5 


13 948. 


2. 9396 


5 


8 



J 9' 1 
270. 49 
270. 49 
270. 50 
270. 50 

270. 48 
270. 48 
270. 48 
270. 48 



abs w 
0. 17523 
. 17474 
. 17417 
. 17418 

.06016 
. 05971 
. 05706 
. 05700 



abs w 
» 0. 17451 
.17411 
. 17419 
. 17430 

. 05694 
. 05680 
. 05709 
. 05709 





percent 


abs w 


1 


f +0. 16 


) 


1 1C- 4-49 


I -0.07 
| -0.02 
I +0.04 

f +0.07 


[ 0. 174231 
| 


i 5-14-50 


J -0.32 

+0.19 

1 +0.19 


\ 0. 05698) 



Days 



137.7 

b (±1.0) 



Sample 2 



10-5-49. . 



4-27-50. 
4-28-50- 



11:10 a. m. 
12:45 p. m. 
2:15 p. m. 

1:48 p. m. 
1:00 p. m. 



6 168. 5 
5 224.0 
5 478.0 


9. 1630 
7. 7486 
8. 1294 


23 
32 
40 


32 

40 

48 


16 514. 9 

17 999. 8 


8. 8470 

9. 5623 


24 
33 


33 
43 



270. 56 
270. 57 
270. 58 

270. 50 
270. 50 



0. 40190 
. 40133 
. 40154 



. 14491 
. 14370 



0. 40183 
. 40139 
. 40173 

. 14425 
. 14373 



4-28-50 



+0.04 
-0.06 
+0.02 

f +0. 18 
I -0.18 



}«• 



14399J 



138. 52 
' (±0. 4) 



Sample 3 



10-8-49 

4-19-50 
4-20-50 



1:24 p. m. 
2:45 p. m. 
4:06 p. m. 

1:49 p. m. 
11:35 p.m. 
2:31 p. m. 



5 134. 4 


12. 2025 





12 


5 207. 5 


12. 3623 


12 


24 


4 576. 2 


10. 8695 


24 


35 


12 692. 3 


11.4817 





11 


10 313. 6 


9. 2898 


11 


20 


10 808. 6 


9. 7279 


20 


30 



270. 58 
270. 61 
270. 63 

270. 51 
270. 52 
270. 52 



0. 643C6 
. 64241 
. 64280 

. 24471 
. 24367 
. 24347 



0. 64323 
. 64278 
. 64335 

. 24358 
. 24365 
. 24360 



4-20-50 



+0.02 
-0.05 
+0.04 

-0.01 

+0.02 

0.00 



0. 24361. 



138. 52 
> (±0. 23) 



Sample 4 



10-6-49.. 



10-7-49. . 



5-12-50. 



12:55 p. m. 
1:35 p. m. 
1:28 p.m. 
2:07 p. m. 
2:57 p. m. 
3:43 p. m. 
4:28 p. m. 
5:14 p. m. 
6:00 p. m. 

11:00 a. m. 
11:54 a. m. 
12:50 p. m. 

1:44 p. m. 

2:35 p. m. 

3:28 p. m. 

10:04 a. m. 
12:25 p. m. 

2:46 p. m. 

5:00 p. m. 

6:55 p. m. 

8:43 p. m. 



2 634. 7 


14. 0020 


28 


42 


2 766. 9 


14. 7066 


42 


57 


2 622. 8 


13. 8785 





14 


2 676. 2 


14. 1607 


14 


28 


2 733. 6 


14. 4582 


28 


42 


2 717. 2 


14. 3735 


42 


57 


2 736. 5 


14. 4643 


57 


71 


2 689. 3 


14. 1840 


71 


86 


2 879. 7 


15.2113 


86 


101 


3 170. 1 


10. 2320 


10 


20 


3 345. 9 


10. 7884 


20 


31 


3 359. 8 


10. 8281 


31 


42 


3 067. 7 


9. 8837 


42 


52 


3 115. 1 


10.0448 


52 


62 


3 178.3 


10. 2384 


62 


72 


8 282. 3 


14. 8070 


3 


18 


8 395.9 


15.0148 


18 


33 


8 534.5 


15. 2193 


33 


48 


7 439. 2 


13. 2613 


48 


61 


6 260. 2 


11. 1775 


61 


72 


6 947. 3 


12. 3852 


72 


85 



270. 84 
270. 91 
270. 72 
270. 78 
270. 84 
270. 91 
270. 97 
271.03 
271. 10 

270. 64 
270. 67 
270. 70 
270. 73 
270. 74 
270. 78 

270. 56 
270. 58 
270. 61 
270. 63 
270. 64 
270. 66 



1. 43937 


1.43994 


1. 43251 


1. 43279 


1. 43249 


1.43306 


1. 43226 


1.42948 


1.43202 


0. 87353 


. 87274 


. 87242 


. 87225 


. 87302 


. 87228 


. 48370 


. 48389 


. 48257 


. 48243 


. 48322 


. 48252 



1. 43307 
1. 43382 
1. 43296 
1.43345 
1. 43338 
1. 43423 
1. 43361 
1. 43105 
1.43382 

0. 87335 

. 87272 
. 87257 
. 87256 
. 87349 
. 87292 

. 48350 
. 48393 

. 48284 



.48391 
. 48339 



N 


/ -0.01 






+0.04 






-0.02 






+0.01 




\ 10- 7-49 


+0.01 
+0.07 
+0.02 
-0.15 
I +0. 04 

f +0.05 
-0.03 


1. 43327\ 

> 


1-14-50 


-0.04 

-0.04 

+0.06 

. +0.00 

f +0.00 
+0.10 


0. 87294\ 

1 


5-12-50 


-0.12 

-0.10 

+0.09 

0.00 


0. 48342J 



138. 39 

b (±0. 14) 



138. 40 
b (±0. 18) 



a Experiment weighted one-third because of short duration of experiment. 

b Figure given for uncertainty based upon authors' estimate of both accidental and systematic errors. 



samples having higher power because it would re- 
quire an increase of only 0.0003 w in the value ob- 
served on May 14 to agree with the other samples. 

It is interesting to note the corresponding effects 
of an uncertainty of 0.0003 w on the other samples 
having higher powers. With samples 2 and 3, this 
uncertainty in the May values would result in half- 
life uncertainties of about 0.21 percent (0.29 day) 
and 0.12 percent (0.17 day), respectively. With 
sample 4, however, the same power uncertainty in 
the values on January 14 results in only 0.04 percent 
(0.06 day) uncertainty in the over-all half -life. 
The uncertainty in heat leak is probably the largest 
error in the measurements with samples 1,2, and 3. 
With sample 4, however, the heat leak uncertainty is 



probably comparable with all other errors. In addi- 
tion, sample 4 was observed at an intermediate date 
(0.7 half -life) in an effort to detect any change in the 
decay constant with time. 

The measurements with the four samples agree 
within ±0.1 percent, or 0.0003 w, whichever is the 
larger. The slightly different half-life values on 
sample 4 for the two different periods are not 
significant because the difference is much less than the 
experimental error. There is no significant evidence, 
therefore, of a change in the decay constant, such as 
would occur either if the polonium was contaminated 
with radioactive impurities having values of half-life 
different from polonium, or if some secondary chemi- 
cal or nuclear reaction produced heat that was not 
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directly proportional to the radioactive energy of the 
polonium sample. In addition, the polonium sam- 
ples used were of a purity that should have precluded 
the possiblity of appreciable contamination by other 
radioactive elements. The polonium had been 
purified both by distillation and by electrodeposi- 
tion on platinum foils from dilute nitric-acid solution 
of polonium nitrate. Of course, the stable lead 
isotope formed from the decay of polonium was pres- 
ent but could not affect the results of this investiga- 
tion. 

The best value of the half-life of polonium was 
estimated from the two values for sample 4, giving 
more weight to the earlier value, where the sample 
had a larger power. On this basis, the authors be- 
lieve that the best value of the half-life derived from 
these measurements is 138.39 days. It is believed 
that the calorimetric uncertainties (based upon both 



accidental and systematic errors) in the experiments 
may result in an error in this value of ±0.1 percent 
(0.14 day). 

The results of these measurements are in agree- 
ment with the value reported by Beamer and Easton 
[1] of 138.3 days ±0.1 percent, determined calori- 
metrically by another method. 
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